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摘 要 : 为 了 探究 陕西 黄土 区 农田 主要 士 壤 养分 特征 及 其 景 
(TN) 有 效 磷 (AP) 和 速效 钾 (AK) 含 量 进行 空 
sx 间 变 异 的 解释 程度 。 结 果 表 明 :SOM TN AP 和 AK 含 量 均值 分 别 为 


数据 ,利用 CIS 与 地 统计 学 方法 对 有 机 质 (SOM) ER 
地 理 探 测 器 模型 探究 18 种 影响 因素 对 养分 空 


区 响 因素 ,基于 研究 区 5096 个 农田 0~20 cem 耕 层 采样 点 
s 间 分 析 ,并 应 用 


14.43 g*kg '.0.92 g- kg 18.21 mg: kg 和 190.28 mg: kg ,呈现 中 等 程度 变异 ;4 种 养分 最 佳 拟 合 模型 均 为 指数 模型 ， 


各 养分 呈现 中 等 程度 空间 相关 性 ,结构 性 与 随机 性 因素 的 共同 作用 导致 了 养分 含量 的 空 
区 域 差异 


关 性 大 小 表现 为 :TN>SOM>AK>AP; 养 分 含量 


明显 ,呈现 


s 间 差异 ;养分 全 局 空间 相 
由 北向 南 逐 渐 递 增 的 趋势 ;年 日 照 时 长 .年 均 气 


温 、 化 肥 用 量 和 地 貌 类 型 等 单 因 子 作 用 对 各 养分 含量 空间 变异 具有 更 强 的 解释 力 ,两 因子 交互 作用 对 养分 的 解释 


力 强 于 其 单 因子 解释 力 。 研 究 表明 , 陕 北 地 区 宜 适 当 增 力 


方面 因素 。 


关键 词 : 土壤 主要 养分 ; 空间 变异 ; 影响 因素 ; 地 理 探测 需 模 型 ; 


农田 土壤 养分 含量 直接 影响 着 农业 生产 效益 
和 粮食 安全 ,是 衡量 耕地 质量 的 重要 指标 H 
中 ,SOM 为 作物 提供 各 类 养分 并 改善 土壤 理化 性 
质 ,N PK 是 影响 作物 生长 发 育 的 关键 养分 ”。 
土壤 属性 、 地 形 地 貌 和 气候 条 件 等 随机 性 因素 m 
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上 肥料 投入 ,关中 地 区 宣 进 行 精耕细作 ,农田 建设 应 考虑 多 


陕西 黄土 区 


研究 。 地 理 探测 器 模型 在 农田 养分 研究 上 的 应 用 
相对 较 少 , 相 比 于 其 他 方法 ,此 方法 无 线性 假设 ,可 
为 土壤 养分 影响 因素 研究 提供 更 合理 的 解释 "”"。 黄 
土 区 作为 黄河 流域 重要 的 生态 单元 ,其 农田 养分 研 
究 备 受 关注 。 尉 芳 等 中 运用 地 统计 方法 对 渭 北 黄 
土 区 农田 土壤 SOM 进行 研究 ,发 现 研 究 区 SOM 总 体 


xd dE ais 3 间 异 质 性 “。 了 解 农田 土 壤 
养分 空间 分 布 .变异 特征 及 其 影响 因素 作用 对 农田 
土壤 健康 可 持续 发 展 及 推动 农业 高 质量 发 展 具有 
重要 意义 
近年 来 ,大 量 学 者 对 农田 土壤 养分 分 布 状况 及 
影响 因素 进行 了 研究。 人 研究 对 象 包 含 SOM NP、K 
FERFI KA Fe, Mn, Cu, Zn 等 微量 元 素 *“ ,人 研 
究 范 围 覆 盖 了 田 块 "果园 "大 小 流域 "各 级 行 
政 区 域 *”“ 等 尺度 。 地 统计 方法 和 GIS 技 术 克 服 了 传 
统统 计 学 解释 养分 空间 变化 上 的 缺陷 ,是 当前 土壤 养 
分 空间 分 析 较 为 有 效 的 方法 之 一 TERI 
尔 逊 相关 性 分 析 外 回归 分 析 中 \、 主 成 分 分 析 "" 方 
差分 析 '* 等 方法 被 广泛 应 用 于 农田 养分 影响 因素 
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上 呈现 南 高 北 低 的 空间 格局 并 揭示 了 其 最 优 预测 
模型 ; 解 文 艳 中 等 对 黄土 高 原 潇 河流 域 农田 土壤 
SOM 进行 了 研究 ,发 现在 该 尺度 下 施肥 等 因素 对 
SOM 含 量具 有 更 强 影响 。 

当前 ,大 多 研究 聚焦 于 中 小 尺度 上 的 SOM ,大 
尺度 下 同时 考虑 区 域 多 类 养分 的 研究 较 少 ,因而 难 
以 为 农田 精准 施肥 提供 更 全 面 的 指导 ,大 尺度 农田 
土壤 养分 分 布 规律 有 待 进一步 揭示 。 此 外 ,大 多 
研究 仅 选择 了 部 分 影响 因素 ,并 未 全 面 考虑 各 类 因 
素 的 影响 。 鉴 于 此 ,本 文 以 陕西 黄土 区 为 研究 对 
象 , 考 虑 多 类 影响 因素 ,人 研究 其 农田 主要 土壤 养分 
特征 及 其 影响 因素 作用 , 拟 解决 以 下 问题 :(1) 了 解 
主要 养分 空间 分 布 格局 ,为 区 域 精准 施肥 提供 理论 
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依据 ;(2) 通过 地 理 探测 器 模型 探究 大 尺度 下 各 类 
因素 与 农田 土壤 养分 的 关系 ,为 陕西 黄土 区 农田 工 
程 建设 及 农业 土壤 养分 管理 提供 参考 。 


1 研究 区 与 研究 方法 


1.1 研究 区 概况 

陕西 黄土 区 位 于 黄土 高 原 中 部 ,地 处 33"42'~ 
39?35'N,106?18'-111?15'E ZR] CE 1) ,高程 范围 为 
207-3154 m ,南北 跨度 大 ,东西 狭窄 ,总 面积 约 13x 
10* km2, 占 全 省 面积 的 3/5, 占 黄土 高 原 面 积 的 1/5。 
区 域 主要 农业 土壤 类 型 包括 陕 北 的 黄 绵 土 和 风沙 
土 . 关 中 平原 的 褐 土 .河谷 滩地 的 潮 士 和 新 积 土 。 
按 地 貌 类 型 ,研究 区 由 北 至 南 划 分 为 :长 城 沿 线 风 
沙滩 区 、 陕 北 黄土 丘陵 沟 蜜 区 、 滑 北 黄 土台 坟 区 和 
关中 平原 区 ; 按 气候 条 件 可 依次 划分 为 :长 城 沿线 
温带 半 干 旱 气 候 区 . 陕 北 高 原 暖 温带 半 干 旱 气 候 区 
和 关中 平原 暖 温带 半 湿 润 气候 区 。 
1.0 样品 采集 及 数据 获取 
1.2.1 养分 指标 选取 与 土壤 样品 采集 SOM TN、 
AP AK 的 丰 缺 程度 是 决定 土壤 肥力 的 重要 因素 , 直 
接 影响 农作物 的 生长 发 育 , 被 广泛 选取 作为 农田 土 
壤 肥 力 人 研究 的 重要 指标 。 此 外 ,4 种 养分 均 属 于 
陕西 省 耕地 质量 评价 中 的 主要 土壤 养分 , 故 本 文 以 
该 4 类 养分 为 研究 对 象 。 结 合 研 究 区 行政 区 位 与 土 
地 利用 现状 ,综合 考虑 地 形 .气候 .土壤 性 质 等 因素 
布设 采样 点 ,于 2020 年 农作物 收获 后 采集 0~20 em 
耕 层 土壤 样品 ,采样 时 用 GPS 记录 点 位 坐标 及 高 
程 ,并 记录 各 样 点 的 土壤 类 型 和 作物 类 型 等 ,本 次 
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共 采 集 5096 份 样品 ,采样 位 置 见 图 1。 样品 采集 及 
养分 检测 工作 由 陕西 省 农业 厅 和 耕地 质量 与 农业 环 
境 保护 工作 站 实施 完成 。 
1.2.2 影响 因子 选择 与 获取 ”综合 以 往 的 研究 成 果 
及 人 研究 区 实际 情况 ,选择 18 种 影响 因子 如 表 1 
所 示 。 

土壤 属性 数据 及 人 为 活动 中 的 灌溉 方式 .作物 
类 型 和 熟 制 数据 来 源 于 陕西 省 农业 厅 耕 地 质量 与 农 
业 环境 保护 工作 站 2020 年 耕地 质量 监测 资料 。 坡 度 、 
坡 向 和 地 形 起 伏 度 通过 空间 地 理 数 据 云 (http:/www- 
Gscloud.cn/) 下 载 的 30 m ASTER GDEM 在 ArcGIS 
10.7 中 提取 ;2020 年 气象 数据 从 中 国 科 学 院 资源 环境 
科学 数据 中 心 (http://www.resdc.cn/) 下 载 ,在 ArcGIS 
10.7 中 提取 采样 点 位 年 日 照 时 长 .年 降雨 量 、 年 平均 
气温 .年 蒸发 量 ; 地 貌 类 型 分 区 通过 陕西 省 农业 农 
村 厅 官 网 (http://nynct.shaanxi.gov.cn/) 发 布 的 关于 农 
田 建设 的 相关 信息 确定 ;区 (县 ) 化 肥 使 用 量 通过 各 
市 区 统计 年 鉴 整理 而 得 。 
1.3 研究 方法 
1.3.1 半 方 差 函 数 分 析 ”地 统计 学 中 的 半 方 差 函 数 
分 析 能 够 揭示 结构 性 因素 与 随机 性 因素 对 空间 养 
分 含量 变化 的 影响 程度 站 。 其 计算 公式 如 下 : 


N(h) , 
y(h)= va; 4402 -Z(*h) — (D 


式 中 : y(h) 表示 半 方 差 函 数值 ;h 表 示 两 采样 点 的 间 
距 ( 或 称 为 步 长 );N(h) 表 示 间 距 为 的 数据 配对 数 ; 
Z(x). Z(x;* h) 分 别 表 示 变 量 Z(x) 在 空间 区 位 x 和 
x+ 天 处 的 实测 值 。 半 变异 函数 的 重要 参数 包括 块 
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图 1 研究 区 区 位 及 采样 点 分 布 
Fig. 1 Distribution of study sites and sampling sites 
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表 1 土壤 养分 影响 因子 
Tab.1 Soil nutrient impact factors 
类 型 变量 影响 因子 单位 /分 类 
土壤 属性 pH X1 一 
土壤 容重 X2 g:cm? 
土壤 类 型 X3 HE KVE 3A He IREE XP E ARAE kE 
成 土 母 质 X4 黄土 母 质 、 风 沙 母 质 、 冲 积 母 质 \ 湖 积 母 质 \ 残 积 母 质 、 坡 积 母 质 \ 红 土 母 质 
土壤 质地 X5 MEDRE EHE TE .中 培土 GEHE E RETE .黏土 
气候 因素 年 日 照 时 长 X6 h 
年 降雨 量 X7 mm 
年 平均 气温 X8 "C 
年 蒸发 量 X9 mm 
地 形 因 素 ”高 程 X10 m 
坡度 X11 
坡 向 X12 阳 坡 (0° ~45°、315° ~360°)、 半 阳 坡 (45° ~ 90°、270° ~ 315°)、 
半 阴 坡 (90° ~ 135°、225° ~270°)、 了 表 坡 (135° ~ 225°)( 以 正 北方 向 为 基准 ) 
地 表 起 伏 度 X13 m 
地 貌 类 型 X14 长 城 沿 线 风 沙滩 区 , 陕 北 黄土 丘陵 沟 蜜 区 、 渭 北 黄土 台 卉 区 、 关 中 平原 区 
人 为 活动 ”化 肥 用 量 X15 kg: hm? 
灌溉 方式 X16 无 灌溉 漫灌 ENE NE TAFE ,喷灌 
作物 类 型 X17 谷类 作物 .果蔬 作物 、 葛 类 作物 、 药 用 作物 \ 豆 类 作物 \ 油 料 作 物 、 其 他 
E X18 年 一 熟 .一 年 两 熟 .一 年 三 熟 .两 年 三 熟 .多 年 生 
DE Co ERE Cte AEA RERANCH Ue 象 背后 各 影响 因素 的 驱动 力 强 弱 "”。 因 子 探测 器 


ERAR FP , 块 金 系数 CV(Co+C) 用 来 反映 变量 
的 空间 相关 性 ”。 

1.3.2 空间 自 相关 分 析 ”为 反映 陕西 黄土 区 土壤 养 
分 在 整体 空间 上 的 聚集 状态 ,采用 全 局 Moran's 7 
指数 进行 计算 分 析 的 。 其 公式 如 下 : 


My>>mile Je E z) 


{LJE 
| 2 


(Xv, s-a) 


P N IRR; 号 为 空间 权重 ;ww “为 1 和 7 处 养 
分 含量 实测 值 ; x 为 养分 均值 。Moran”s 1 指数 值 域 
为 [-1,1j, 用 以 表示 养分 空间 分 布 的 相关 性 ,通常 
采用 Z 得 分 检验 指数 显著 性 。 

1.3.3 空间 聚 类 分 析 通过 局 部 Moran’ s 7 指数 能 
够 反映 区 域 局 部 空间 自 相 关 性 从 而 进行 养分 空间 
聚 类 分 析 和 空间 异常 值 识别 六 ,公式 如 下 : 


N 
N(x; 一 DAA 
JA 
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1.3.4 地 理 探测 器 模型 ”地 理 探测 器 模型 以 统计 学 
手段 对 数据 空间 分 异性 进行 探测 ,并 有 效 地 描述 现 


I- (2) 


N 
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i=l 
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用 于 探测 影响 因子 X 对 因 变 量 Y 的 解释 程度 (解释 
力 ) ,度量 指 标 用 g 统 计量 表示 ,gq 的 值 域 为 [0,1], 值 


越 大 ,因子 解释 力 越 强 ,其 公式 如 下 : 
YN; 
1-5. -1-30 (4) 
WP RR A EUR, 2…, L); 07, 为 影 


WATES h ESSERI EL AR 8 B 77 28 5 o^ 为 因 变量 在 
全 范围 的 方差 ;WN 为 所 有 因子 的 数量 ;WN; 为 第 h 层 时 
因子 数量 ;SS 下 为 因子 层 内 方差 之 和 ;SST 为 总 体 方 
差 之 和 ;4 为 因子 对 因 变 量 了 的 影响 力 。 

交互 探测 器 描述 两 因子 共同 作用 对 因 变 量 的 
解释 程度 。 通 过 装 加 两 因子 分 层 重新 计算 g 值 ,将 
两 因子 共同 作用 的 g 值 与 单 央 子 g 值 对 比 ,判断 交互 
作用 强 弱 。 根 据 对 比 , 将 交互 作用 分 为 5 类 ,如 表 2 
所 示 。 
1.4 数据 处 理 

对 5096 份 采样 数据 采用 3 倍 标准 差 法 与 邻近 
点 数据 比较 法 吻 除 各 类 土壤 养分 数据 异常 值 , 别 除 
异常 值 后 SOM ,TN .AP、AK 样 本 数 分 别 为 5079 份 、 
5084 份 4995 份 .5020 份 。 通 过 SPSS 27.0 进 行 统计 
性 描述 ,通过 GS+9.0 进行 半 变 异 函 数 分 析 , 在 Arc- 
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R2 双 因 子 交互 作用 类 型 据 的 偏 斜 程度 , 峰 度 系数 表示 数据 的 陡峭 程度 , 偏 
Tab.2 Type of two-factor interactions 度 越 趋 于 0, 峰 度 越 趋 于 3, 数 据 越 近似 服从 正 态 分 
判断 依据 交互 作用 f ^" ,对 AP 与 AK 进 行 对 数 转化 后 ,4 种 养分 含量 均 
qGOGUG)»q (X )q QG) 非 线 性 增强 近似 服从 正 态 分 布 。 
q 17Td 2J]2q 1 2J]2max| 143 2 ü 时 z DE 4 
OONAN ATA 22 土壤 养分 空间 变异 及 分 布 特 全 
q i d0445/—q i +q 2 ESZ iol B ` 
X o É kiA y 壤 养 分 空间 变异 4. ` 
asiste IO srai 2.2.1 基于 半 方 差 模型 的 土壤 养分 空间 变异 特征 通 
min[g CX.) ,gq(X,)] 减弱 过 GS+9.0 对 研究 区 养分 分 布 状况 进 和 了 半 方 差分 析 ， 
q (XN) min[g (Xi) ,q(6)] 非 线性 持 搞 结果 如 表 4 所 示 。 分 析 结 果 表 明 ,4 种 养分 的 最 佳 拟 


GIS 10.7 中 进行 土壤 养分 空间 自 相 关 分 析 、 养 分 空 
间 插 值 和 空间 聚 类 分 析 ,通过 有 R 语 言 构建 地 理 探 测 
器 模型 探讨 各 因素 对 土壤 养分 含量 的 影响 。 


2 结果 与 分 析 


土壤 养分 统计 特征 

土壤 养分 统计 特征 分 析 ( 表 3) 表 明 ,SOM TN、 
AP 和 AK 含 量 均值 分 别 为 14.43 g .kg ,0.92 g.kg ， 
18.21 mg .kg ,190.28 mg*kg'。 与 全 国 二 次 土壤 普 
查 养分 分 级 标准 进行 对 比 ,SOM 与 TN 处 于 区 等 ( 较 
低 ) 水 平 ,AP 与 AK 分 别 为 亚 等 (中 等 ) 工 等 ( 较 高 ) 
水 平 。4 种 养分 变异 系数 分 别 为 :42.27% ,42.3996 、 
79.90% 和 56.88%。 根 据 Nielsen(1985)Xl] Z) ERME", 
C. 0.1 为 弱 变 异 ,0.1 < Cs 1 为 中 等 程度 变异 强度 ， 
C>1 为 强 变异 ,各 项 土壤 养分 指标 均 属 于 中 等 变异 
水 平 ,AP 和 AK 相 对 于 SOM 和 TN 变异 强度 较 高 ,说 
明 两 者 存在 更 明显 的 区 域 差异 。 在 后 续 地 统计 分 


2.1 


合 模型 均 为 指数 模型 。SOM .TN AP 、AK 的 尼 均 趋 
近 于 1 , 拟 合 效果 优良 。 各 养分 块 金 值 Co 均 大 于 0， 
这 是 由 于 采样 误差 .实验 误差 及 最 小 取样 间距 内 的 
人 为 随机 因素 所 导致 的 空间 变异 ”。 块 金 系 数 Cw 
(Co+O) 表 示 随 机 变量 导致 的 变异 在 总 变异 中 所 占 的 
比例 ” ,4 种 养分 的 块 金 系数 位 于 0.25~0.75 之 间 , 表 
明 养 分 含量 呈现 中 等 程度 相关 性 ,变异 性 同时 受 土 

壤 属 性 、 气 候 条 件 、 地 形 条 件 等 结构 性 和 人 为 活动 
类 随机 因素 的 影响 。AK 块 金 系数 为 0.68, 相 对 偏 
大 ,说 明 其 受 人 为 活动 影响 相 较 其 他 养分 更 大 。4 
种 养分 分 别 在 40.8 km、58.9 km、92.4 km .67.3 km 的 


范围 内 具有 相关 性 。 
222 土壤 养分 空间 自 相 关 分 析 Moran'sI 8X 


大 于 0 表示 养分 空间 分 布 为 正 相 关 ,在 空间 上 呈 聚 
集 趋 势 ; 小 于 0 表示 养分 空间 分 布 为 负 相 关 ,在 空间 
上 旦 分 散 趋势 ;等 于 0 表示 养分 含量 之 间 不 相关 , 空 
间 上 趋 于 随机 分 布 。 在 0.01 置信 水 平 下 ,IZI> 1.96、 
2.58 分 别 表 示 养 分 含量 呈现 显著 自 相 关 性 、 极 显著 


析 中 需要 对 数据 进 4 TIER, ,否则 可 和 g 导 致 变异 
函数 分 析 和 空间 插值 的 偏差 中 。 偏 度 系数 表示 数 


表 3 陕西 黄土 区 土壤 养 


AHR, WAKEN Moran’ s 1 指数 及 标准 化 Z 
值 如 表 4 所 示 , 各 养分 全 局 Moran' s 1 指数 值 均 大 于 


量 统计 特征 


Tab.3 Statistical characteristics of soil nutrient content in the Loess Plateau of Shaanxi Province 


指标 样本 数 最 大 值 最 小 值 平均 值 标准 差 变异 系数 /% 分 布 类 型 WE MEE 
SOM/Cg*kg') 5079 32.97 1.26 14.43 6.10 42.27 正 态 分 布 -0.10 2.52 
TN/Cg* kg) 5084 2.10 0.04 0.92 0.39 42.39 正 态 分 布 -0.18 2.50 
AP/(ng-*kg*) 4995 71.28 0.04 18.21 14.55 79.90 对 数 正 态 -0.42 3.50 
AK/(mg:kg'') 5020 550.00 18.00 190.28 108.24 56.88 对 数 正 态 0.04 2.56 


表 4 半 方 差 消 数 模型 参数 及 全 局 Moran s 7 指数 


Tab.4 Semi-variance function model parameters and the global Moran’ s 7 index 


HER 论 模 型 REC 基 台 值 C(,t+C” 块 金 系数 Cy(CotC) ” 变 程 4/km  UUOEARGC 残 差 平方 和 RSS 全 局 Moran s7 Z 
SOM H n 0.1941 0.3892 0.4987 40.8 0.979 3.79x10^ 0.58 86.31 
TN 指数 模型 0.0134 0.0231 0.5931 58.9 0.995 3.31x10* 0.61 90.90 
AP 指数 模型 0.3530 0.7070 0.4993 92.4 0.994 5.56x10™ 0.31 44.79 
AK 指数 模型 0.1573 0.2299 0.6842 67.3 0.964 2.46x10* 0.36 53.40 
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0, 且 IZI 值 大 于 2.58, 表 明 4 种 土壤 养分 均 呈 现 极 显 


风沙 滩 区 呈 焉 级 标准 ,含量 中 等 ;位 于 关中 平原 地 


著 正 相关 , 即 某 一 小 区 域 土 壤 养分 含量 值 高 低 会 受 
到 周转 同一 养分 含量 值 高 低 的 正 影响 。 养 分 全 局 


区 的 西安 宝鸡 渭南 三 市 及 其 周边 区 域 2 种 养分 含 


Moran”s 1 指数 值 表现 为 :TN(0.61) >SOM(0.58) > 
AK(0.36) > AP(0.31) , TN SOM 空间 聚集 现象 相对 
于 AP 和 AK 更 明显 。 

2.2.3 土壤 养分 含量 空间 分 布 特征 ”根据 半 方 差 模 
型 拟 合 所 得 结果 参数 ,采用 普通 克 里 金 法 进行 空间 
插值 ,通过 自然 间断 点 法 将 土壤 养分 含量 分 为 5 个 
级 别 ( 图 2)。4 种 养分 含量 总 体 区 域 差异 明显 ,呈现 
南 高 北 低 的 特征 。SOM 和 TN 的 空间 分 布 特征 相 
似 ,这 是 因为 SOM 的 积累 与 分 解 对 氮 元 素 的 存储 和 
转化 有 着 直接 影响 ,因而 对 土壤 氮 素 含量 有 着 主导 
作用 2。 榆 林地 区 及 延安 北部 两 者 含量 相对 匮乏 ， 
多 处 于 JV 、V 级 标准 ,延安 南部 和 榆林 西 侧 的 局 部 


(a) SOM 


图 例 
有 机 质 含量 /(g-kg 1!) 
| «9.59 

[33 9.59~12.91 
E] 12.92~16.41 
ES 16.42-20.00 
EH 220.00 


m 


(c) AP 


>z 


图 例 

有 效 磷 含量 /mg-kg-) 
ES <9.03 

[3 9.03-14.35 
E] 14.36-20.74 
国 国 20.75-29.62 
EN >29.62 


VOTI 


EREET I. IRRE, S ETE ,并 向 南 呈 现 递 
增 趋 势 。AP 含 量 V 、V 级 标准 区 域 面积 较 大 ,其 中 ， 
V 级 标准 多 位 于 延安 及 榆林 东 侧 , T REEL 
分 布 于 西安 市 南 侧 边 缘 。AK 含量 由 西北 向 东南 呈 
现 递增 趋势 , 低 值 分 布 在 榆林 和 延安 西 侧 , 高 值 分 
布 于 东南 平原 区 域 。 

224 土壤 养分 空间 聚 类 分 析 ”基于 普通 克 里 金 插 
值 结果 ,以 (区 ) 县 为 评价 单元 ,将 (区 ) 县 内 插值 结 
果 的 平均 值 空间 连接 到 评价 单元 ,并 进行 局 部 Mo- 
ran"s 了 指数 计算 ,形成 LISA 聚 类 图 (图 3)。 空 间 聚 
类 分 析 将 研究 区 划分 为 :不 显著 .高 -高 型 .高 - 低 
型 . 低 - 高 型 \. 低 - 低 型 5 类 区 域 ,分 别 表示 区 域 之 间 
高 低 值 无 显著 相关 性 ,高 值 聚 集 .自身 高 值 但 周围 


(b) TN 


图 例 
PAR E/(gkg') 
E «0.60 

[3 0.60-0.79 
E] 0.80-1.01 
ERI 1.02-1.24 
EH >1.24 


>z 


(d) AK 


图 例 
速效 钾 含 量 /mgkg-) 
BE «117.69 

[3 117.69-166.08 
E] 166.09-215.98 
ER 215.99-270.43 f 
EH 27043 


m 


Ai Aer 


2 土壤 养分 空间 分 布 


Fig.2 Spatial distribution of soil nutrients 
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(a) SOM 


(c) AP 


图 3 陕西 黄 


区 农田 了 


-Z 


图 例 

M 不 显著 
EB 高 -高 
EB 高 - 低 
EB 低 - 高 
EB 低 - 低 


0 


100 km 


(d) AK 


À 


JANT 


EIRIAN LISA 聚 类 分 布 


Fig.3 LISA cluster distribution of farmland soil nutrients in the Loess Plateau of Shaanxi Province 


低 值 自身 低 值 但 周围 高 值 . 低 值 聚 集 。 该 方法 常 
以 LISA 聚 类 图 反映 研究 区 局 部 土壤 养分 之 间 的 关 
联 性 的 。 结 果 表 明 ,SOM 与 TN 聚集 特征 较为 相似 ， 
低 值 聚 类 区 均 分 布 于 榆林 地 区 及 延安 北部 ,高 值 聚 
类 区 分 布 于 西安 市 及 其 周边 ,TN 高 值 聚 类 区 也 包括 
了 宝鸡 市 大 部 分 区 域 。AP 低 值 聚 类 分 布 于 研究 区 
中 西部 地 区 ,高 值 聚 类 区 多 分 布 于 西安 市 及 其 周 
边 ,在 咸阳 市 北部 有 小 部 分 高 - 低 值 区 。AK 低 值 聚 
类 区 分 布 于 榆林 及 延安 西北 范围 ,高 值 聚 集 区 多 分 
布 于 西安 .成 阳 渭南 三 市 。 
2.3 土壤 养分 影响 因素 分 析 
2.3.1 单 因 子 作用 探测 “通过 因子 探测 器 计算 18 
影响 因子 对 SOM TN AP 和 AK 空 间 分 布 的 q 统 
计量 以 表示 其 解释 力 , 除 坡 向 因素 对 4 种 养分 含量 
的 影响 力 不 显 著 (P>0.01) 外 ,其 他 因子 均 达 到 显著 


水 平 (P<0.01)。 因 子 探测 结果 (图 4) 表 明 , 各 影响 
因子 对 区 域 养分 空间 分 异 的 解释 力 (g 值 ) 存 在 明显 
差异 。 总 体 来 看 ,各 因子 对 SOM 和 TN 的 解释 力 显 
著 大 于 AP 和 AK, 这 说 明 各 因素 在 农田 SOM 与 TN 
循环 过 程 中 参与 度 相 对 于 AP 和 AK 更 高 。 就 单一 
影响 因素 而 言 , 对 SOM 解释 程度 最 高 的 5 类 因子 依 
次 为 :年 日 照 时 长 (0.44) .地 貌 类 型 (0.42) ,年 气温 
(0.35) 化 肥 用 量 (0.34) 作物 熟 制 (0.30) ,对 TN 解 
释 程 度 最 高 的 5 类 因子 依次 为 :年 日 照 时 长 (0.46)、 
化 肥 用 量 (0.37) 地貌 类 型 (0.36) ,年 气温 (0.32) 4E 
蒸发 (0.26); 对 AP 解释 程度 最 高 的 5 类 因素 分 别 
为 :年 气温 (0.17) ,高 程 (0.17) .地貌 类 型 (0.16) .化 
肥 用 量 (0.14) .年 日 照 时 长 (0.14) ;对 AK 解 释 程 度 
最 高 的 5 类 因子 分 别 为 :年 气温 (0.25) .年 日 照 时 长 
(0.24) , & ££ (0.22) 化 肥 用 量 (0.21) 、 地 貌 类 型 
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(0.20)。 对 4 类 养分 解释 位 于 前 5 的 因素 均 包 含 了 
年 日 照 时 长 .年 均 温度 、 地 貌 类 型 化肥 用 量 。 
通过 交互 探测 右 探 究 各 影响 
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图 4 陕西 黄土 区 农田 养分 影响 因子 g 值 


Fig. 4 q-value of nutrient influence factor in the Loess Plateau of Shaanxi Province 


2.3.10 交互 作用 探测 


因素 对 土壤 养分 含量 属于 独立 作用 还 是 相互 作用 ， 
如 果 是 相互 作用 ,是 增强 还 是 减弱 。 列 出 各 养分 前 
5 的 交互 作用 因子 如 表 5 所 示 。 

根据 表 5 np, 


因子 的 交互 作用 对 SOM 5j TN zs 
间 分 异 的 解释 力 更 强 , 位 于 前 5 的 g 值 均 大 于 0.45; 
而 对 于 AP 和 AK, 因 子 交互 作用 的 影响 稍 弱 , 分 别 位 
于 0.22 与 0.30 左 右 。 对 SOM 解 释 最 大 的 交互 作用 


表 5 不 同 土壤 养分 因子 交互 作用 探测 


Tab.5 Interaction detection of different soil nutrient factors 
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因子 为 pH 与 年 均 日 照 时 长 (0.496), 大 于 pH(0.18) 与 
年 均 日 照 时 长 (0.44) 单 因子 作用 ,说 明 两 因子 交互 
作用 对 养分 的 解释 相 比 于 单 因子 作用 的 解释 力 显 
著 提 升 。 对 TN 解释 最 大 的 交互 作用 因子 与 有 机 质 
TRIR] Cg IEN 0.535) ,这 一 方面 是 因为 两 因素 对 TN 含 
量 的 直接 影响 , 另 一 方面 是 其 作用 于 SOM ,间接 的 
影响 了 TN 的 空间 差异 。AP 与 AK 受 年 均 温 度 与 化 
肥 用 量 的 交互 影响 最 大 ,g 值 分 别 为 0.229 5 0.303, 
说 明 人 为 因素 与 自然 因素 的 交互 对 两 养分 有 更 强 
的 解释 力 。 交 互 因 子 作 用 热 图 (图 5) 表 明 ,各 类 因 


" SOM TN AP AK 
交互 作用 RONDE NEUES SMS s 

q 交互 类 型 q 交互 类 型 q 交互 类 型 q 交互 类 型 

交互 作用 1 X1NX6 双 因 子 X1NX6 WAF X8NX15 双 因 子 X8NX15 双 因 子 
0.496 曾 强 0.535 兽 强 0.229 增强 0.303 曾 强 

交互 作用 2 X14NX15 双 因子 X6NX15 双 因 子 X10NX14 双 因 子 X10mX15 双 因 子 
0.492 兽 强 0.512 Eh 0.227 增强 0.299 兽 强 

交互 作用 3 X6NX18 HAF X2NX6 双 因 子 XImX8 双 因 子 X6NX10 双 因 子 
0.487 兽 强 0.506 Bh 0.226 增强 0.298 兽 强 

交互 作用 4 XImX12 双 因 子 X5NX6 双 因 子 XImX10 双 因 子 X7mX8 双 因 子 
0.479 兽 强 0.505 Bh 0.225 增强 0.295 DES 

交互 作用 5 XImX15 双 因 子 X6NX9 双 因 子 X15NX16 双 因 子 X6NX8 双 因 子 
0.478 兽 强 0.498 兽 强 0.222 增强 0.294 兽 强 
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图 5 土壤 养分 因子 交互 作用 


Fig. 5 Soil nutrient factor interaction diagram 


子 对 陕西 黄土 区 农田 养分 的 交互 作用 解释 力 均 高 
于 单 因子 作用 ,表现 为 双 因 子 增强 与 单 因子 非 线性 
增强 。 


3 讨论 


空间 格局 及 变异 特征 方面 ,统计 分 析 表 明 , 陕 
西 黄土 区 农田 土壤 养分 整体 偏 低 ,其 原因 与 区 域 复 
杂 的 土壤 、 气 候 、 地 形 等 自然 条 件 和 粗放 的 耕作 方 
A ,管理 制度 等 人 为 因素 密切 相关 ”。AP 的 变异 系 
数 为 79.90% ,显著 高 于 其 他 养分 ,这 是 因为 磷 元 素 
易 被 土壤 固定 上 且 更 容易 受 人 为 施肥 灌溉 等 因素 影 
啊 ”。 从 理论 来 看 ,陕西 黄土 区 南北 地 形 地 貌 . 气 
候 差 异 明 显 ,土壤 养分 应 受 结构 性 因素 影响 程度 更 


高 ,但 半 方 差 函 数 模型 的 结果 表明 ,4 种 养分 均 呈 现 
中 等 程度 相关 性 ,养分 含量 受 结构 作用 与 随机 作用 
共同 影响 ,导致 这 一 差异 的 主要 原因 可 能 是 由 于 长 
期 的 环境 保护 与 农业 栽培 管理 使 土壤 养分 发 生 了 
一 定 的 变化 ,因而 随机 因素 的 影响 逐渐 加 强 ”。 养 
分 的 空间 格局 分 布 与 聚 类 状况 表明 ,榆林 与 延安 北 
部 4 类 养分 缺乏 ,这 可 能 是 由 于 该 地 区 多 为 风沙 土 
与 黄 绵 土 ,土壤 本 身 贫 交 , 且 环境 条 件 使 得 养分 流 
失 较 为 严重 ”。 秸秆 还 田 及 施用 有 机 肥 不 仅 可 以 
为 农田 土壤 养分 提供 丰富 的 有 机 质 及 毛 磷 钾 , 还 可 
以 改善 土壤 结构 ,适合 在 当地 推广 ""。 西 安 、 宝 鸡 、 
渭南 一 带 ,地 处 关中 平原 , 测 土 . 褐 土 等 土壤 养分 会 
量 丰 富 , 且 地 势 相 对 平坦 ,土壤 保水 保 肥 能 力 强 ， 
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而 多 为 养分 富 集 区 ”。 
影响 因素 方面 ,土壤 属性 通过 决定 保水 保 肥 、 


上 壤 主 要 养分 特征 及 影响 因素 1915 


在 农田 建设 过 程 中 ,应 注意 影响 作物 产量 的 多 方面 
因素 ,并 根据 不 同 地 区 土壤 属性 ` 地 形 条 件 .气候 条 


通气 和 耕作 难 易 程 度 影响 土壤 养分 含量 ”。 土 壤 
质地 相 比 于 其 他 土壤 属性 解释 力 更 高 ,这 是 因为 质 
地 较 轻 的 土壤 通 透 性 较 大 ,更易 发 生 SOM 的 矿 化 和 
其 他 养分 的 流失 ,不 利于 养分 积累 ,而 质地 较 重 的 
土壤 往往 竹 粒 含量 高 ,更 容易 吸附 和 固定 SOM 等 养 
分 “。 气 候 条 件 往往 可 以 决定 养分 空间 大 格局 ,日 
照 时 长 .降水 量 温度 、 蒸 发 等 通过 影响 作物 光合 
作用 、 区 域 水 热平衡 微生物 作用 进而 影响 养分 含 
量 "”。 其 中 ,日 照 时 长 作用 最 为 明显 ,人 研究 区 内 养 
分 含量 总 体 上 随 日 照 时 长 的 降低 而 增高 。 地 形 地 
貌 主要 通过 影响 水 分 与 太阳 辐射 的 空间 再 分 配 从 
而 影响 养分 含量 ,其 中 ,地 貌 类 型 与 高 程 在 一 定 区 
域内 控制 着 气候 土壤, 水文 等 因素 的 地 域 性 差异 ， 
是 制约 农田 养分 的 重要 因素 ""。 人 为 活动 方面 ,化 
肥 用 量 作为 耕地 养分 的 直接 来 源 , 对 养分 含量 空间 
差异 具有 更 强 的 解释 力 ,而 耕作 模式 .管理 制度 等 
因素 也 间接 影响 着 养分 含量 。 交 互 探测 表明 两 因 
子 交 互 作 用 对 养分 的 影响 均 强 于 单 因子 影响 ,说 明 
了 养分 含量 受 各 因素 的 共同 作用 ,这 与 王 琦 等 ”、 
尉 芳 等 ”在 该 区 域 的 研究 结果 一 致 。 


4 结论 


(1) 根据 全 国 第 二 次 土壤 普查 养分 分 级 标准 ， 
陕西 黄土 区 整体 SOM 与 TN 处 于 较 低 水 平 ,AP 与 
AK 分 别 为 中 等 水 平和 较 高 水 平 , 各 土壤 养分 指标 
均 属于 中 等 变异 水 平 。 

(2) 4 种 养分 最 佳 拟 合 模型 均 为 指数 模型 ,各 养 
分 呈现 中 等 程度 空间 相关 性 ,结构 性 与 随机 性 因素 
的 共同 作用 导致 了 研究 区 农田 土壤 养分 含量 的 空 
间 差 异 ; 全 局 Moran' sx 指数 表明 ,4 种 养分 均 呈 现 
极 显著 正 相 关 ; 克 里 金 插值 表明 ,SOM , TN, AP 和 
AK 含 量 区 域 差异 明显 , 呈 南 高 北 低 的 特征 ;空间 聚 
类 分 析 表 明 , 低 值 聚集 区 多 分 布 于 延安 榆林 一 带 ， 
高 值 聚 集 区 多 分 布 于 西安 .宝鸡 .渭南 一 带 。 建 议 
陕 北 地 区 适当 增加 肥料 投入 ,多 采取 保水 保 肥 措 
Ji ,关中 地 区 进行 精耕细作 ,提高 粮食 产量 和 养分 
利用 效率 。 

(3) 年 日 照 时 长 .年 均 气 温 、 化 肥 用 量 和 地 貌 类 
型 对 各 养分 含量 空间 变异 特性 解释 力 更 强 , 两 因子 
交互 作用 对 养分 的 解释 力 均 强 于 其 单 因 子 作 用 。 


件 等 因素 因地制宜 的 进行 人 为 活动 。 
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Abstract: To explore the main soil nutrient characteristics and influencing factors of farmland in the Loess 
Plateau of Shaanxi, spatial analysis was conducted on the content of SOM, TN, AP, and AK using GIS and 
geostatistics methods based on the data from 5096 farmland sampling points with a depth of 0-20 cm in the study 
area. A geographic detector model was used to explore the degree to which 18 influencing factors explained the 
spatial variation of nutrients. The results showed that the average contents of SOM, TN, AP, and AK were 14.43 
g:kg ', 0.92 g- kg ', 18.21 mg: kg ', and 190.28 mg-kg ', respectively, showing moderate variation. The four best- 
fit models for nutrients were all exponential models, with moderate spatial correlation among each nutrient. The 
combined effect of structural and random factors causes spatial differences in the nutrient content. The global 
spatial correlation of nutrients is TN > SOM > AK > AP. The regional differences in nutrient content are 
significant, indicating a gradually increasing trend from north to south. The single- factor effects of annual 
sunshine duration, annual average temperature, fertilizer use, and geomorphic type have stronger explanatory 
power for spatial variation in nutrient content but lesser than that of the interaction between the two factors. 
Research has shown that it is necessary to increase fertilizer input in the northern Shaanxi region, perform 
intensive cultivation in the Guanzhong region, and consider multiple factors in farmland construction. 

Keywords: main nutrients in soil; spatial variation; influencing factors; geological detector model; the Loess 
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